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The observation of photons with energies above 10 18 eV would open a new window in 
cosmic-ray research, with possible impact on astrophysics, particle physics, cosmology 
and fundamental physics. Current and planned air shower experiments, particularly the 
Pierre Auger Observatory, offer an unprecedented opportunity to search for such photons 
and to complement efforts of multi-messenger observations of the universe. We summarize 
motivation, achievements, and prospects of the search for ultra-high energy photons. 
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1. Introduction 

Photons are the main messenger particles for exploring the universe. Over the last 
decades, the wavelength range of photon observation was dramatically expanded and 
stretches now from radiowaves to high-energy gamma rays. The current maximum 
energy of photons observed is ^10 14 evP 

In this short review, we focus on the search for photons of energies above 10 18 eV 
up to the highest energy (few times 10 20 eV) by measuring particle cascades (air 
showers) initiated in the atmosphere of the Earth. An observation of such ultra- 
high energy (UHE) photons, several orders of magnitude in excess of currently 
observed photon energies, would open a new window in cosmic-ray research with 
significant impacts on related research fields. Giant air shower experiments, most 
prominently the Pierre Auger Observatory^ 1 , are unique tools to explore this photon 
energy range. The expected sensitivity of the Auger Observatory after installation 
of a large northern array reaches the photon fluxes predicted in more conservative 
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scenarios of cosmic-ray origin, making an observation of UHE photons feasible. 

UHE photons a re also inherently lin ked to other messenger particles such as 
charged cosmic rays^^EEl an( j neutrinoJ^l The search for UHE photons com- 
plements current experimental efforts towards multi-messenger observations of the 
universe. 

The structure of the paper follows the possible "life" of UHE photons. Produc- 
tion scenarios and propagation of photons as well as their predicted fluxes at Earth 
are shortly described in Section [2] Geomagnetic and atmospheric cascading and 
the challenge of identifying photon showers are discussed in Section [3l Current ob- 
servational results and future prospects are summarized in Section [4] Throughout 
Sections [2j-|4l conceptual issues important for the data interpretation are pointed 
out. The possible impact of photon searches and detections on different research 
fields is briefly outlined in Section [5l Conclusions are given in Section [6] 

2. UHE Photons As Cosmic Rays 

Production. Though at different levels, most models predict UHE photons, mainly 
from the decay of neutral pions produced previously by a "primary process" , 

primary process -> tt° (+tt ± ) + ... -> 7 UHE (+i> UHB ) + ... (1) 

In non-acceleration models^, the primary process i s given by the decay or annihila- 
tion o f primordial relics such as topological defect er^^ l (TD) or super heavy dark 
matter™™ 

(SHDM). From considerations of QCD fragment ation! 16 | 17 | 18 | 19 | 
copious photons are then expected to be producedPSEHZSl j n the Z-burst 
scenario22124125| (ZB), photons are generated via the resonant production of Z 
bosons by UHE neutrinos annihilating on the relic neutrino background. 

In more "conventional" cosmic-ray models, nuclear primaries are accelerated at 
suitable astrophysical sites to ultra-high en ergy. U HE photons can be produced 
during propagation by the GZK-type proces a 1 1 of resonant photo-pion produc- 
tion of UHE nucleons with the cosmic microwave background. The energies of these 
"GZK photons" are typically a factor ^10 below the primary nucleon energy. An 
enhanced production of photo ns > 10 18 eV can occur from nuclear primaries passing 

IOQI 

near the galactic center region. 1 ^ 

Propagation. UHE photons can initiate electromagnetic cascades by interacting 
with background radiation fields, 

7lJHE ^background * ^ > > 7gcV-TcV •" (2) 

Significant uncertainties exist for the low-frequency (fe w MH z) radio background 
and (for e ± propagation) extragalactic magnetic fieldsPQEU Typical energy loss 
lengths assumed for UHE photons range between 7-15 Mpc at 10 19 eV and 5- 
30 Mpc at 10 20 eV. In SHDM models, the relic particles are clustered as cold dark 
matter in our Galaxy, and UHE photons, as their decay products, would be observed 
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Fig. 1. Energy loss length of photons for 
interactions with infrared (IR), cosmic mi- 
crowave (CMB) and universal radio (URB) 
backgrounds. Uncertainties exist for the 
URB (see text) and IR background. For 
comparison, curves for protons and iron nu- 
clei are add ed. Values are compiled from 
Ref s [9l29p0l, 



at Earth with little processing. In TD and ZB models, UHE photons are injected 
at larger distance from the Earth, and part of these photons can cascade to lower 
energy. The electromagnetic cascade stops at GeV-TeV energies where the universe 
becomes increasingly transparent for photons, see Figure [1] 

Flux predictions. In addition to the mentioned uncertainties of UHE photon 
propagation, there are free parameters within the theoretical source models such as 
density and lifetime of primordial relics in non-acceleration models, injection spec- 
trum in acceleration models, source distribution, etc. To calculate the (range of) 
photon fluxes predicted by a model, the measured cosmic-ray spectrum is usually 
taken as a constraint. The question arises which spectrum to use. Firstly, the ab- 
solute UHE cosmic-ray flux differs between the experiments by a factor ^2. This 
problem may be reduced by regarding the fraction of photons in the cosmic-ray flux 
rather than absolute r ates. S econdly, it is unclear so far whether a suppression of 
the flux ("GZK cutoff'^SEll) a b 0V e E Gm - 6 x 10 19 eV exists (as indicated by 
HiRes datcPJ) or not (as indicated by AGASA data^). The sha pe of t he assumed 
energy spectrum can affect also the predicted fraction of photons BH33] 

Non-acceleration models are usually calculated assuming a spectrum without 
flux suppression. These models aim at explaining the highest-energy end of the 
spectrum onlyP2l Yoy energies up to ~-Egzk , a "conventional" component of nuclear 
primaries must be assumed. Then, a common prediction of the models is a dominant 
photon component at ~10 20 eV, see Figure [2] A reduction of the predicted photon 
fraction may formally be achieved by introducing a flux suppression above Eqzk 
which increases the relative importance of the "conventional" nuclear component. 
However, at the same time this reduces the relevance of non-acceleration models to 
explain any observed events at all. 

For acceleration models, in turn, it seems more natural to assume a spectrum 
with flux suppression as no nearby astrophysical sources could be identified by now. 
The predicted photon fluxes are typically relatively small. Fractions of order ^0.1% 
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Fig. 2. Fraction of photons in the integral cosmic-ray flux as a function of the threshold en- 
ergy. The predictions of non- acceleration models (dashed lines: ZB, SHDM, TD from ED, SHDM' 
fromESl) follow a spectrum without flux suppression above Egzk- Experimental l imit s fSection l4.ll 
all limits at 95% c.l. ) are r eported from Auge r hy brid obse rvations (Auger hybridj331, Haverah Park 
measurements (HP^EESl AGASA data f AfP ZI. fA2| SH], Yakutsk data (Yp, preliminary j22| and a 
combination of AGASA and Yakutsk data (AY) 10 41 I To illustrate the minimum number of events 
required to place a (95% c.l.) limit to a certain photon fraction, horizontal dotted lines are shown 
with the number of events assigned. An estimate (Section 14.21 of the sensitivity of the southern 
Auger Observatory for two years of operation is given for the two cases that the real spectrum 
(I) does have, (II) does not have a cutoff. The data from the southern Auger Observatory allow a 
stringent test of non-acceleration scenarios even if the real spectrum had a cutoff. 



were obtained in scenarios assuming nucleon sources j I I with a considerable range 
of uncertainty (see Figure [2]). Larger photon fractions may emerge for specific as- 
sumptions on source features, particularly when trying to reproduce a spectrum 
without flux suppression!^ In case of primary nuclei, the rate of GZK produced 
photons may be suppressed due to the higher total energy of the nucleus required. 
Further investigations of such scenarios and of the relation between the photon flux 
and model parameters are desirable. 

A caveat seems in place. The reconstruction of the energy spectrum from shower 
observations itself requires, to some extent, the fraction of photons as an input (Sec- 
tion [JO]). Usually, no contribution from primary photons is assumed in the recon- 
struction. For model predictions, a self-consistent comparison to data is required. 

3. UHE Photons as Air Shower Primaries 

At the Earth, photons initiate almost purely electromagnetic showers via pair pro- 
duction and bremsstrahlung. Additional processes are important at highest energy. 

3.1. Specific high-energy processes 

Preshower effect. Contrary to nuclear primaries, ^10 20 eV photons can con- 
vert in the geomagnetic field to an pair which then emits synchrotron 
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Fig. 3. Sky maps of photon energies at which the conversion probability IS -Pconv = 50% for the 
southern (left) and planned northern (right) Auger sites I^SI Contour lines are given with a stepsize 
of Alg(£/eV) = 0.1, minimum and maximum energies are assigned. Azimuthal directions are 
labeled ("E" for East etc.). Zenith angles are given as concentric circles in 10° steps (zenith in the 
center). With increasing energy, the sky fraction grows where the preshower process is important. 
At Auger North, preshowering starts at factor ~2 smaller energy compared to Auger South due 
to the stronger local magnetic field. Close to the pointing direction of the local magnetic field 
(indicated by a red dot), the highest energies are required for conversion to occur. Accordingly, 
the sky patterns are shifted in local coordinates between the sites. Despite the larger magnetic 
field at Auger North, higher energies are needed to reach P C onv = 50% close to the local field 
direction. This is connected to the field lines being less curved with altitude at Auger North. The 
differences between the site characteristics can be exploited in photon searches 1^21 

photonsP B 4 | 45 | 46 | 47|48 B51 Instead of a single UHE photon, a bunch of lower- 
energy electromagnetic particles, called "preshower", enters the atmosphere with 
important consequences for the shower development. The local differential proba- 
bility of photon conversion as well as the probability distribution of synchrotron 
photons emitted by the electrons, depend on the parameter 



where E is the energy of the parent particle (photon or electron), m the electron 
mass, B c a constant, and B± is the local magnetic field component transverse to the 
direction of the particle's motion. Since B± is involved, preshower characteristics 
depend strongly on the specific trajectory through the magnetosphere and, thus, on 
the arrival direction and the experimental site (see Figure [3]). The probability P C onv 
of a photon to convert in the Earth's magnetic field results from an integration 
along the particle trajectory; for P conv — > 1, photons would almost surely undergo 
geomagnetic cascading. Non-negligible probabilities (P CO nv ~ 10%) are usually ob- 
tained if values \ > 0.5 are reached along the trajectory, corresponding to photon 
energies above 2-4 x 10 19 eV depending on the site. 

The energy partition in the produced e ± pair is usually symmetri c (rat io of 
energies <2), with the probability of asymmetric share increasing with ^ HSUE The 
spectrum of synchrotron photons becomes harder for larger values of \- Secondary 
photons of sufficiently high energy can convert again. 



X = 



E B± 
mc 2 B c 



B c ~ 4.414 x 10" G 



(3) 
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A typical preshower from a 10 20 eV photon starts at ~1000 km altitude and 
enters the atmosphere (~100 km altitude) with one or a few e ± pairs around 10 18 eV 
and a large number (~500) of photons. The photon energies extend over several 
decades, with a few photons around 10 19 eV. Significant fluctuations around these 
averages may occur, particularly when the conversion takes place at low altitudes. 
The spread of preshower particles in transverse distance and arrival time is well 
belo w cu rrent detector resolutions; the subsequent air shower is observed as one 
event 

To simulate the preshower effect, Monte Carlo codes have been developed. The 
PRESHOWER^ program is available as a standalone tool and is also linked to the 
CORSIKA^and CONEX 54 shower codes. An independently devel ope d preshower 
program (MaGICEp^) is available within the shower code AIRESfSSI. A realistic 
model of the geomagnetic field such as the IGRP^l model, is an essential ingredient 
of precise calculations. As an example, a simple dipole model fails to reproduce the 
factor ^2 difference in local field strength at the southern and northern Auger sites. 

LPM effect. In a medium, the Bethe-Heitler cross-sectional for pair produc- 
tion by photons (ctbh ~ 0.51 b in aiJ^Sj) can be reduced due to destructive in- 
terference fr om sev eral scattering centers. This so-called LPM effect (Land au and 
PomeranchuliSQISE MigdaP^J) is confirmed by experiments, see 

Refe.lMlMI f or 

re- 
views. With 

E~,E LPM „ m 2 c 3 aX Q . X 

Elpm = — -7—^ ~ ( 7 - 7 TeV/cm) x — , (4) 



E e {E 1 -E e y urLVL Anhp v ' ' ' p 

the reduced cross-section olpm can for k<1 be approximated^^ by ulpm = 
&bh\/~k oc (pE^)~2^ with photon energy E 1 , electron energy E e , radiation length 
Xq ~ 37 g cm -2 in air, density p of the medium, and electron mass m. It follows 
from Eq. ((4|) that the reduction is largest for conversion to a symmetric electron pair 
(E e w E 1 /2). The cross-section for producing a highly asymmetric pair (E K /E 7 — > 
or 1) changes only slightly. In a similar way, also bremsstrahlung is suppressed ESESI 
Numerical examples (see e.g. Ref.^j) are Slpm - 2.8 x 10 17 eV at 300 m a.s.l. and 
~ 10 19 eV in the upper atmosphere. 

The LPM effect delays the development of an air shower initiated by a single 
UHE photon as those processes are suppressed that degrade the energy carried by 
individual high-energy particles. Fluctuations can be very large due to a positive 
correlation of the reduction of olpm since <jlpm{X2) < <Jlpm ( Ai ) for depths 
A 2 >Ai. The LPM effect is accounted for in AIRES, CORSIKA, and CONEX. 

Photonuclear interactions. Photon- initiated cascades are almost purely elec- 
tromagnetic ones. Production of muon pairs is suppressed by {m e /m^) 2 . The cross- 
section for photonuclear interactions which mainly transfer energy to secondary 
hadrons (and these subsequently to muons), is expected to be ~10 mb at 10 19 eV 
and thus more than two orders of magnitude below <7bh (which, however, can be 
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Fig. 4. Observed profile of the 
3.2 X 10 20 eV Fly's Eye evenJZS 
(data po ints ) compared to sim- 
ulations assuming primary 
photons when neglecting both 
prcshower and LPM effect (dot- 
ted line), accounting for LPM ef- 
fect only (dotted-dashed line), and 
accounting for both effects (solid 
line). Data points are correlated in 
X. Photons differ by ~1.5o" from 
the data which correspo nds to a 
~13% chance probabilityESES 



reduced by the LPM effect). 

One important consequence is that our lack of knowledge of hadron dynamics at 
high-energy — a major limitation for conclusions about the nuclear composition — 
is of much smaller impact when calculating photon shower s. F or instance, depths 
of shower maxima differ by 30—40 g cm -2 between SIBYLLrS3 and QGSJET 68 69 

for UHE protons (see Section [3~2l Figure [5]) but by less than ~ 5 g cm -2 for UHE 

1041 

photons^. For photon searches, data can be compared to photon simulations only, 
without subtracting a background from nuclear primaries. 

As another consequence, certain photon shower observables, in particular the 
number of secondary muons, are sensitive to UHE extrapolations of the photonu- 
clear cross-section. Different e xtrapo lations of the photon-proton cross-section exist 
in the literature, see e.g. RefsP^^ for compilations. For UHE photon searches, it 
is important to know whether the cross-section could gro ssly ex ceed an extrapola- 
tion such as provided by the Particle Data Group (PDG) E2H3J smce larger values 
mean photon showers were actually more similar to those of nuclear primaries; an 
absence of UHE photons may then erroneously be concluded from data. Using a 
dipole formalism and taking unitarity constraints into account, values of the UHE 
cross-section that exceed the PDG extrapolation by more than ~80% (less at smaller 
energy) were found to be disfavored!^ It should be noted that for shower simula- 
tions, finally the photon- ai r cr oss-section is needed, which requires an account for 
nuclear effects (see e.g. RefP^). From simulations with modified cross-sections, un- 
certainties for photon showers of ~10 g cm -2 in X max (depth of shower maximum) 
and ~15% in iV M (number of muo ns at ground) were estimated!^ Assuming an 
extrapolation such as from Refill which exceeds the PDG fit by a factor ~10 at 
10 20 eV, a reduction of A max by ~100 g cm -2 is possible for unconverted photons. 
Already with a small sample of observed photon showers, such a scenario can be 
distinguished from more "standard" ones. 
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Fig. 5. Average depth of shower maxi- 
mum <X max > versus energy simulated 
for primary photons, protons and iron 
nuclei. The impact of the LPM and 
preshower effects on <X 

max is visi- 
ble. The splitting of the photon line 
at ~3xl0 19 eV indicates that <X max > 
above these energies depends also on 
the specific trajectory through the geo- 
magnetic field. For nuclear primaries, 
calculations for different hadronic inter- 
action models are dis played (SIBYLL 
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See Refill for references to the experi- 
mental data. 



3.2. Features of photon-induced showers 

Shower profiles calculated with the preshower and LPM effects switched on/off are 
shown in Figure [4] for the conditions of the 3.2 x I0 20 Fly's Eye eventl^. Both 
effects compete with each other: the LPM effect increases the average X max and its 
event-by-event fluctuations. The preshower effect, in turn, reduces both quantities. 

Calculations of average X max values for different primaries are shown in Figure[5] 
The large "elongation rate" (slope dX max /d lg E) for photons leads to X max values 
well above those of nuclear primaries already at I0 17 — 10 18 eV. Changes in the 
elongation rate at 10 19 — fO 20 eV are due to the LPM and preshower effects. 

In addition to X max , showers from photon and nuclear primaries differ in due 
to the small photonuclear cross-section. To study the impact of shower fluctuations 
on photon identification, a scatter plot of X max versus is given in Figure [6] for 
different primaries simulated at 3 x 10 19 eV. This is an idealized plot as no detector 
effects are accounted for. It can be seen that shower fluctuations by themselves do 
not impose major limitations to identify photons. Misidentification rates of <10~ 3 
(<10 -4 ) using X max (Nfj) seem possible for protons and are even smaller for nu- 
clei (numbers may depend on the hadron generator used for calculating nuclear 
primaries). In turn, an unconverted primary photon may hardly be distinguished 
from a nuclear primary if in one of the first interactions, a photonuclear interac- 
tion occurs (probability of order Contrary to X max , does not differ much 
between converted and unconverte d photo of photon showers 
can be found, for instance, in Refs! a l 0U l 01 l 0/ l 'J l 04 l 00 l 0l:) l I. 

There are other observables, usually closely related to X max or iV M or a combi- 
nation of both, that could be used to distinguish photons from nuclear primaries. 
These observables, some of which are specific for a certain detector type and experi- 
mental configuration, i ncl ude the curvature of the shower fron t and the steepness of 
the lateral distribution^ f ground particles, the risetime^of the detector signal 
at a certain core distance as well as the signal amplitude by itself. 
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Fig. 6. X max versus JV M for dif- 
ferent primaries at 3 X 10 19 eV. 
Simulations were performed for 
photons (1000 events), protons 
(10000 events), and iron nuclei 
(100 events) with PRESHOWER, 
CONEX and QGSJET 01. No de- 
tector effects are included. 
The lines illustrate possible cuts 
to achieve small misidentification 
rates for photon identification (see 
text). 



3.3. Detector response to photon-induced showers 

The differences between showers from photons and nuclear primaries can affect the 
detector acceptance to and the energy reconstruction of UHE photons. As the effects 
depend on the detector type, analysis approaches appropriate for the experimental 
apparatus can be developed to minimize uncertainties from the corrections. 

Suppose a two-component composition of photons and protons with fluxes ^ 1 (E) 
and & P (E) at energy E. The ratio of photons to protons, one measure of the photon 
fraction, is then f lp (E) — $>~ f (E)/$ p (E). Similarly, integral fluxes or fractions, 
additional primaries, or the photon fraction to the total flux could be considered. 

Regarding the detector acceptances Ai (E) with i — 7 or p, an acceptance ratio 
of e(E) — Ary(E)/A p (E) < 1 affects both the flux and fraction of registered photons. 
Such a bias can occur for near- vertical photon showers reaching ground before being 
fully developed. For instance, the aperture of the HiRes-I telescope was found to 
be ~40% reduced for photons at Eqzk^^ In a recent analysis of Auger hybrid 
data, a minimum zenith angle was introduced to reduce the acceptance bias to 
photons. 1 ^ Experimental conclusions about the photon flux require knowledge of 
A 1 {E) (usually well controlled for arrays using appropriate selection cuts). For 
conclusions about the photon fraction, knowledge of the relative acceptance e(E) 
(possibly easier to calculate for fluorescence telescopes) may be sufficient. The latter 
approach has been applied in the analysis of Auger hybrid data!^31 

Reconstructing an energy E rec from observables that, on average, differ between 
photon and proton showers of same energy E and assuming proton primaries only, 
leads to a misreconstruction of photon energies. Selecting events according to E Tec 
then introduces a shift in energy scales between the primaries: events with the same 
E rcc refer to protons of energy E ICC = E, but to photons of energy E TCC /g = E/g 
with g^l. Hence, the photon flux and fraction actually entering the data sample 
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are $> 7 (E/g) ^ &-y(E) and & 1 (E/g)/$ p (E) ^ f lv {E). In general, a correction for 
this shift in energy scales is required. 

For the fluorescence technique using the integrated calorimetric shower energy 
as an observable, g ~ 1.1, i.e. the energy scales nearly match and g > 1. The 
small difference in energy scales is due to the missing energy correction, which is 
larger for nuclear primaries compared to the almost purely electromagnetic photon 
showers 1^21 A conservative upper limit to the integral flux or fraction of photons 
above E may then be possible without correction or assumptions on $ 7 (_E) sinc e for 
g ~ 1.1 => E/g <E^ J E/g $ 7 (E')dE' > J E %{E')dE' . This was used in Ref. ^ 

For arrays when reconstructing the energy from a signal amplitude S(r) mea- 
sured at ground at core distance r, the factor g depends in general on the arrival 
direction of the event (slant depth to ground changes with zenith, preshower proba- 
bility depends on zenith and azimuth) , its energy (dependence of A max with energy, 
see Figure 03), and dete ctor ch aracteristics (e.g. sensitivity to shower muons). Values 
of g 0.5 are possible! 11 Additionally large shower fluctuations may occur (for 
a proposal to reduce the effect of shower fluctuations when reconstructing photon 
energies, see Refill). The (for arrays larger) uncertainty when relating different 
energy scales may be avoided by analysing the photon flux rather than the fraction, 
as for the flux, only the photon energy scale is required (and A 1 {E) is usually well 
controlled for arrays, see above). 

Within the flux level experimentally allowed for photons, the different detector 
response to photons introduces a systematic uncertainty when reconstructing an 
energy spectrum. An example of how the high-energy end of the spectrum from 
first Auger data can change when simply assuming a constant factor two between 
photon and proton energies is show n in RefP. Detector responses to deep showers 
in general are discussed in Refp^l 

4. Experimental Searches for UHE Photons 
4.1. Status 

No photon detection has been reported so far. Upper limits to photons come from 
different experiments, see also Figure [2] 

Comparing rates of near-vertical showers to inclined ones recorded with the 
Haverah Park water detectors, upper limits (95% c.l.) of 48% above 10 19 eV (52 
events) and 50% above 4x 10 19 eV (10 events) were deducedPSESlFrom simulations, 
ground signals were found to be stronger suppressed for photons than for nuclear 
primaries at large zenith angles. The absence of this suppression in the data gives 
constrains on the photon contribution. Contrary to other approaches, a feature of 
the data sample as a whole is analysed here instead of observables in single events. 

With muon counters in the AGASA array, the muon density at 1000 m core 
distance was measured for part of the events!^ A simulation study of the ground 
signal S(600) for different primaries showed a possible underestimation of photon 
energies (e.g. about 30%, 50%, 20% at 10 19 eV, 10 19 5 eV, >10 20 eV). Assuming a 
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mixture of primary photons and protons, limits (95% c.l.) to the photon fraction 
were estimated to be 28% above 10 19 eV (102 events) and 67% above 3.2 x 10 19 eV 
(14 events) The data were compared to an overall simulated distribution, i.e. not 
event-by-event. There seem to be some "photon-like" events (not commented on in 
the paper): in about five events, the measured p^ is considerably (factor ~4 or more) 
below typical values expected for nuclear primaries; in additional four events, the 
(at least two) muon detectors registered no signal; one event of E Tec ~ 7.5 x 10 19 eV 
(energy scale for nuclear primaries) was observed with a value p^ about a factor 
~2.5 below an average expected from a data fit. 

A comparison to individual events is most interesting at highest energy where 
every second event could be a photon (cf. Figure [2]). Observed event features can 
be compared to high-statistics photon simulations, and the chance probability of 
photons to generate such an event can be determined. This way, for the 3.2 x 10 20 eV 
Fly's Eye evenlP^l ( see Figured]) with X^f x ~ 815 ± 60 g cm -2 , photon shower 
profiles (with J* Q, ~ 937 ± 26 g cm -2 ) were found to differ by ~ 1.5cx from the 
observed profile ] ' ' While a pho ton or igin can not be excluded, profiles from 
nuclear primaries fit the data better jSlZSI 

The method of comparing individual events was also applied in RefPSl to the six 
highest energy AG AS A events with observed p M . The energy scale for photons and 
the uncertainti es o f reconstructed energy and p^ were taken according to AGASA 
findings in RefP3. For all six events, the muon densities from photons are (factors 
2—7) below the observed ones. Photon predictions would be even more discrepant 
to data if smaller energies (see below) were assumed. A statistical approach for 
deriving limits from a small number of events was developed which accounts for 
shower features changing for each event and for non-Gaussian shower fluctuations. 
An upper limit of 67% (95% c.l.) above 1.25 x 10 20 eV was derived. Though not ruling 
out non-acceleration models, these scenarios appear to have problems to consistently 
explain the AGASA data (no flux suppression and no photon dominance). 

A method for comparing measured and simulated shower obs ervable s such as 
S(600) instead of using the reconstructed energy was given in 

Refs ]40141J For indi _ 

vidual events, the corresponding photon energies are found by requiring the simu- 
lated S(600) to fit the observed ones and properly weighting the events. This way, 
shower and detector fluctuations can more d irectly be accounted for. Detailed in- 
formation about the detector is needed32ED 

Applied to AGASA data, using response functions for S(600) and information 
on the muon detector response from RefP, photon energies of the six highest- 
energy events were found in RefP to be (up to factor ~2) below those expected 
from 

RefP. 

There appear to be differences between the findings of Refs!^ and ^3 
in the expected S(600) for fixed energy which may be connected to the use of 
different simulation codes. Combining six AGASA^ZI and four Yakuts]**^ events 
with E ICC > 8 x 10 19 eV (energy scale for nuclear primaries), a 36% limit above 
10 20 eV was derived (95% c.l.)- Th ere may be a sensitivity to the choice of E rec . 
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Reducing E lec , a more "photon- like" AGASA event (see above) enters the sample. 
For a discussion about combining data from these two experiments which app ear 

IOOI 

to have systematic differences in t he reco nstructed flux spectra, see e.g. Ref!^. 

Using the same method of Ref^SEl anc j investigating Yakutsk data only, pre- 
liminary (not yet published) limits of 12% above 2 x 10 19 eV (50 events) and 22% 
above 4 x 10 19 eV (95% c.l.) were obtained!^ In two events no muon signal was 
registered, another event was classified as muon-poor. 

All these limits refer to the fraction of photons and come from ground arrays. The 
fluorescence technique provides important cross-checks due to smaller uncertainties 
from different energy scales of the primaries (Section 13. 3|) . Using A max from the 
direct observation of the shower profile with fluorescence telescopes in hybrid events 
(i.e. registered by ground and fluorescence detectors) a limit of 16% (95% c.l.) above 
10 19 eV (29 events) was obtained from first data taken at the Auger ObservatoryPSI 
An approach was developed such that the limit does not rely on assumptions about 
input spectra or composition. A ^2xl0 20 eV event registered with X°^^ ~ 821 ± 
36 g cm~ 2 differs by ^3<r from the photon hypothesis (A^ ax ~ 948 ± 27 g cm -2 ). 

Steps were taken to analyse the HiRes stereo data set of about 50 events above 
4 x 10 19 eVJ^ No results were reported so far. 

In conclusion, current data from ground and fluorescence detectors do not in- 
dicate a large (above 10 19 eV) or dominant (above 10 20 eV) photon contribution. 
Non-acceleration models are constrained by existing limits though not ruled out as 
an explanation of UHE cosmic rays (see also discussion in Refl^j). 

4.2. Prospects 

The 16% limit from the Auger hybrid data is based on 29 events. The sensitivity 
can be improved by a factor ~3 or more with data accumulated until 2008/2009. 
Also energies below 10 19 eV can be reached with the hybrid technique. 

For an overall estimate of future sensitivities it is instructive to consider the 
lowest theoretical fraction F™ m that can be excluded with n events. With a being 
the confidence level of rejection, F™ m is 

F 7 mil » = 1 - (1 - a) 1/n « 3/n for a = 0.95, n > 1 . (5) 

A limit derived from data can exceed F™ m (i.e. be weaker) depending on the 
discrimination power of the specific observables (possibly limited by shower fluctua- 
tions and detector characteristics), efficiency corrections, and the actual number of 
photons in the data sample. For the example of the Auger hybrid limit based on 29 
events, F™ in (29) - 10%. The final limit of 16%, which includes a - 20% efficiency 
correction, corresponds to the theoretical limit for 17 events (F™ m (17) ~ 16%). 
Thus, compared to the 29 events in the data set, the 16% limit corresponds to an 
effective number of events that is smaller by a factor r ~ 29/17 ~ 1.7. 

To estimate the sensitivity to photons of the large Auger array, we assume 
a factor of r~2.4, i.e. larger (less discriminative events and/or larger efficiency 
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Fig. 7. Fraction of photons in the integral cosmic-ray flux as a function of the threshold energy 
(see also Fig. [2} . The predictions^^ (labeled "GZK photons" ) for primary nucleon sources follow 
a spectrum with flux suppression above Eqzk ■ For reference, the current limit from Auger hybrid 
data, obtained with 29 events, is shown. The estimated sensitivity (see text) of the array of the 
southern Auger Observatory is shown for two and twenty years of operation for a spectrum with 
flux suppression. Also shown is the estimated sensitivity using hybrid data, which are less numerous 
due to the ~10% duty cycle but reach to smaller energy. The uncertainty of the sensitivity estimates 
is indicated at the lower left corner. At around 10 19 eV (dashed lines), additional threshold effects 
may become increasingly important for the array. If complemented by an extended northern array, 
a sensitivity level of below 0.1% can be reached within a few years of full operation. 

correction) than the factor reached in the first Auger hybrid limit. Corresponding 
numbers f rom other event-by-ev ent p hoton searches using AGASA and Yakutsk data 
are r~2.2pSl~2.lp21 and ~1.15^. More detailed studies, particularly of threshold 
effects, are needed. However, the uncertainty of this estimate (factor ~2) seems 
well below that of current theoretical predictions of UHE photon fluxes. Using for 
n(E > Eq) the Auger energy spectrum^ gives the sensitivity estimates F(E > Eq) 
shown in Figures [2] and [7| by calculating F = F" uu (n') ~ 3r/n where n' = n/r 
denotes the effective number of eventsE] The curve F(Eq) indicates the upper limit 
that can be set at 95% c.l. if there were no photons in the data. The curve represents 
also (since F(Eq) ~ 3r/n) the ratio of 3r ~ 7 to n events: if it is possible to identify 
7 photons out of n(Eo) events, F(Eq) gives the observed photon fraction. 

A crucial test of non-acceleration models is possible with the southern Auger 
array even if the flux is suppressed above Eqzk (Figure [2|) . A large exposure is 
required to reach sensitivities of ~ 0.1% or below, i.e. to possibly observe GZK 
photons expected from UHE nucleon sources (Figure [7|). Such sensitivity levels 
may be reached on acceptable time scales if the southern Auger Observatory is 
complemented by an expanded (e.g. factor 3—4 larger) northern site. 

a With the integral number of events n(E > Eo), from Eq. {5} it follows _F™ m oc S c f 1 for a 
power law spectrum of differential index /3 (/3 ~ 2.84 from Ref . 1221) , i.e. a straight line of slope 
(3 — 1 in double-logarithmic scales. Zenith angles were restricted to 30—60° for the estimate of the 
sensitivity. 
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The Telescope Arraj^H f scintillators, under construction in Utah, will cover 
a ~760 km 2 area (~25% of the southern Auger array) overviewe d by fluorescence 
telescopes. Space missions such as the proposed EUSO^^ and OWl]i^ projects may 
offer a significant increase of the experimental exposure to cosmic rays at highest 
energy. 



5. Possible Impact of UHE Photon Searches 

Photons, as the gauge bosons of the electromagnetic force, at such enormous energy 
can be regarded as unique messengers and probes of extreme and, possibly, new 
physics. Implications are related to the production of photons, their propagation, 
and interactions at the Earth. Many aspects of the following, incomplete list of 
possible impacts of UHE photon searches and connections to other research subjects 
require more study. 

Large UHE photon fluxes are a smoking gun for current non-acceleration models. 
Stringent photon limits give p aram eter constraints such as a lower limit on the 
lifetime of relic SHDM particlesPSl 

Findings on photons are needed to reduce corresponding systematics in other 
air shower studies, such as for the energy spectrum (Section 13. 3p or when tryi ng to 
constrain interaction parameters such as the proton-air cross-sectior J^Q^ I ^ Q^ from 
showers. 

UHE photons may be helpful for diagnostics of sources accelerating nuclear 
primaries, as the photon fluxes from UHE hadron interactions are expected to be 
connected with source features such as type of primary, injection spectrum, possible 
beam dump at the source, or source distribution (see also RefPB). 

UHE photons point back to the location of their production. Possibly, the arrival 
directions of photons may correlate better with the source direction than those of 
charged primaries. There may be an enhanced UHE photon flux from the galactic 

IOQI 

center region depending on the spectra of nuclear primaries. 1 ^ 21 In certain SHDM 
scenarios, an enhanced flux of ^10 18 eV p hotons from the galactic center is possible 
without a higher-energy counterpart ! 

Propagation features of UHE photons are sensitive to the MHz radio 
back groun d.^ The photon flux at Earth is also sensitive to extragalactic magnetic 



fields.™ 

Already a small sample of photon-induced showers may provide relatively clean 
probes of aspects of QED and QCD at ultra-high energy via the preshower process 
and photonuclear interactions (Section 13. ljl . 

Certain implications will be most powerful in combination with other results 
from shower observation (spectrum, anisotropy, nuclear composition) and, probably, 
with results using other messenger particles (UHE neutrinos, lower energy photons). 
Parallel to UHE photons, UHE neutrinos are usually produced, see Eq. (TTJ). Due 
to the much longer mean free path of UHE neutrinos, searches for UHE photons 
and UHE neutrinos complement each other, with UHE photons (neutrinos) testing 
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more local (distant) production sites. The disappearance of UHE photons during 
propagation is accompanied by an appearance of GeV-TeV photons, cf. Eq. ([2]) and 
Figure [1] (see also 

Ref.™ 

and references therein). Similar to (but independent of) 
UHE neutrinos, these GeV-TeV photons allow a test of more distant production 
sites. The close relation between different messengers is reflected by the fact that 
constraints on non-accele ration models also come from UHE neutrino and GeV 
photon data, see e.g. RefsP^^H The full information could be exploited by multi- 
messenger observations. For instance, for known sources of UHE cosmic rays, fluxes 
of UHE photons much below the expected level could indicate ce rtain new physics. 

There are several connections to Lorentz invariance violationPIIini The pro- 
duction of GZK photons can be affected as well as interactions of photons during 
propagation and when initiating a cascade at the Earth. Particularly, photon conver- 
sion (interaction with background fields or preshower process) may be suppressed. 

It is interesting to check whether UHE photon propagation could be affected by 
the presence of axions or scalar bosons. Formal requirements for photon co nversion 
regarding photon energy and magnetic field strength, appear to be fulfilled P^DJ^ 
Photon conversions to non-electromagnetic channels may differ from the standard 
QED process of Eq. |(5J) due to an absence of an electromagnetic sub-cascade. 

UHE photon propag ation can be modified in c ertain mo dels of brane worldspH 
quantum gravity theoryP^HHl or spacetime foaiJn3IIII3 p or instance, depending 
on fundamental length scales significant scattering of photons on structures (defects) 
in spacetime foam can occur. In turn, constraints may be d erived wh en actually 
observing UHE photons, even with one gold-plated event only!H41U5| 

6. Conclusions 

UHE photons are tracers of highest-energy processes and new physics. Showers 
initiated by such photons can well be distinguished from those by nuclear primaries. 
When complemented by a large northern site, the Pierre Auger Observatory is 
expected to be sensitive to photon fractions of 0.1% and below and can realistically 
aim at photon observations. The search for UHE photons contributes to multi- 
messenger observations of the universe. A key characteristic for the progress in 
astrophysics is to expand to photon wavelengths beyond the optical; the observation 
of UHE photons would bring this to the highest end. 
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